Coherent optical control of spin-spin interaction in doped semiconductors 
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We provide a theory of laser-induced interaction between spins localized by impurity centers in 
a semiconductor host. By solving exactly the problem of two localized spins interacting with one 
itinerant exciton, an analytical expression for the induced spin-spin interaction is given as a function 
of the spin separation, laser energy, and intensity. We apply the theory to shallow neutral donors 
(Si) and deep rare-earth magnetic impurities (Yb) in III-V semiconductors. When the photon 
energy approaches a resonance related to excitons bound to the impurities, the coupling between 
the localized spins increases, and may change from ferromagnetic to anti-ferromagnetic. This light- 
controlled spin interaction provides a mechanism for the quantum control of spins in semiconductors 
for quantum information processing; it suggests the realization of spin systems whose magnetic 
properties can be controlled by changing the strength and the sign of the spin-spin interaction. 
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I. INTRODUCTION 

The possibility of creating semiconductor systems that 
can work simultaneously as electronic, photonic and mag- 
netic devices has boosted recently the research on light- 
spin interaction in semiconductors. These efforts extend 
both in the direction of an optical control of macroscopic 
magnetic properties and toward the quantum control of 
single spinsii Examples belonging to the first class in- 
clude the investigation on light-induced paramagnetic 
to ferromagnetic transitions in magnetic semiconductors, 
using coherent^ or incoherent^ processes, and the at- 
tempts to control the macroscopic spin polarization of 
carriers using polarized light More recently, schemes for 
optical control at the level of few spins have been pro- 
posed. Bao et ali^ have demonstrated quantum spin en- 
tanglement of few donors and magnetic impurities. The 
strong potentialities of this optical quantum control of 
spins for quantum information processing using quan- 
tum dotsSil and impurities^ have been emphasized, and 
recent advances in the optical control at the nanoscale 
in semiconductor nanostructuresSiifi are particularly en- 
couraging in this perspective. 

It was pointed out in Ref.^that itinerant excitons, i.e. 
optical excitations free to move in the host material that 
embeds the localized spins, can induce an effective spin- 
spin interaction between localized spins. This mechanism 
has been dubbed Optical RKKY (ORKKY), in analogy 
to the mechanism in theory of magnetism^ where elec- 
trons are involved. In the coherent optical case virtual 
excitons are created, and the ORKKY coupling is ob- 
tained from a second order perturbation theory in the 
exchange coupling between the itinerant exciton and the 
localized spin. The ORKKY result predicts that the cou- 
pling between the localized spins is always ferromagnetic, 
independently of the sign of the coupling with the exci- 
tons. In this paper we show that higher order terms in 
the exciton-impurity coupling can modify the strength 
and sign of the interaction, and affect its dependence 
as a function of the spin separation. The calculation of 
the spin-spin interaction can be reduced to a spin de- 
pendent scattering problem that can be solved including 



exactly all the multiple scattering terms between the two 
localized spins. We follow here an approach similar to 
that used to calculate multiple scattering effects of i:^ 
mesons by deuterons^ The higher order terms in the 
exciton-impurity coupling describe bound states which 
affect strongly the optically induced spin-spin interac- 
tion. In particular, a controlled anti-ferromagnetic (AF) 
coupling can be realized when the laser energy is tuned in 
the bonding-antibonding gap for the exciton localized by 
two impurities. This laser controlled switching of sign of 
the spin-spin interaction opens to new directions in the 
investigation of competing interactions in spin systems. 

The extension of the RKKY to insulators was orig- 
inally proposed by Bloembergen and RowlandJ^ They 
predicted the exponential decay of the spin-spin inter- 
action with a characteristic length n ~ 'hj ^ 2mEg de- 
pending on the energy gap Eg and mass of the virtual 
electron hole pairs across the gap. In the optically in- 
duced RKKY, the energy gap is effectively reduced by 
the laser field which increases the effective length to 
K = h/ ^j2m{Eg — fiujp), with hiop being the energy 
of the laser. Also, the density of electrons in the oc- 
cupied bands in the insulator is replaced in the optical 
case by the density of photons in the field. The innova- 
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FIG. 1: Scheme of the light induced spin-spin interaction in 
the case of two shallow donors. 
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tive strength of the optically induced case resides in the 
control potentialities since both the intensity and the fre- 
quency of the laser can be controlled in an experiment. 
There are intrinsic difhculties to go beyond second or- 
der perturbation theory in the case of metals.^'' These 
difficulties are not present in the optical coherent case 
since there is no Fermi sea of electrons. The presence of 
a Fermi sea simultaneously with the laser would produce 
light-induced Kondo effects ji^ which we do not consider 
here. 

The paper is organized as follows: In section |n] we de- 
rive the expression for the effective Hamiltonian of two 
localized spins in the presence of a light field, and we re- 
late it to the spin-dependent T matrix operator of a two- 
center scattering problem. We study first in Sec. Illll the 
scattering of one exciton with one center. By generalizing 
a result from scattering theorjti^ to the spin dependent 
case, we show in section lTvl how the T matrix operator for 
the exciton scattering on two centers can be expressed in 
terms of the T operator for the one center scattering. We 
also study in this section the effects of the polarization of 
the light and we show that a circularly polarized field will 
induce an additional term representing a magnetic field. 
The theory is applied in section to two systems: shal- 
low donors, and deep rare earth magnetic impurities. We 
discuss implications for quantum computing implemen- 
tations and for optical control of macroscopic magnetic 
properties in Sec lVIl 

II. EFFECTIVE HAMILTONIAN FOR 
LOCALIZED SPINS COUPLED BY THE LIGHT 

A schematic view of the effect in the case of two shal- 
low donors is given in Fig. ^ We are not interested in 
calculating the optical properties of the whole system, 
but we want to consider the effect of a coherent field on 
the dynamics of the two non-interacting localized spins 
and . The light creates virtual/real excitons in the 
semiconductor host and affects the localized spin states. 
Due to their coupling to the light, the localized spins can 
therefore interact with each other. We want to study 
the behavior of the two localized spins in the coherent 
optical regime. This implies that the laser is always off 
resonance with respect to the free exciton band to avoid 
strong energy absorption. We therefore consider only sin- 
gle exciton processes in the presence of a monochromatic 
laser field. The system of two localized spins coupled to 
one itinerant exciton is described by the Hamiltonian 

Hxs — Hq + Hi , (1) 

where. Ho describes a free exciton of mass M with dis- 
persion efc = eo + h^k^/2M. The term Hi describing the 
spin dynamics can be written in the form 

Hi = — ^ Jk,k'{s^ ■ Sa'a 
kk' aa' (3 

+ e-^('='-^-)^s^-So-„)6l.,,,.6fc„;3 (2) 



where the two localized spins | are described by and 
. V is the volume, and s is the electronic spin of the 
itinerant exciton. creates an exciton with center 

of mass momentum fc, electron spin a, and hole spin f3. 
R is the separation between the two impurities. Jk,k' 
is the exciton-spin exchange interaction. The strength 
and the sign of this term depend strongly on the na- 
ture of the localized spin. The sign, for instance, is de- 
termined by the competition between the ferromagnetic 
potential exchange and the antiferromagnetic kinetic ex- 
change which is due to the hybridization of the itinerant 
exciton state with the localized state^i We will keep for 
the moment a general approach independent of the na- 
ture of the Jk,k' J and we will discuss two specific exam- 
ples in section A spin-independent term correspond- 
ing to a direct Coulomb interaction between the exciton 
and the impurity is also present. This term is small for 
shallow impurities, where kinetic exchange effects dom- 
inate, but becomes important for deep impurities. We 
will include this term in the case of rare earth impurities 
discussed in Sec. IV Bl and wc disregard it in the general 
discussion since it only introduces spin-independent en- 
ergy shifts. We assume that the Is excitons dominate the 
light induced effect, as discussed in Ref.0 Moreover, we 
focus on systems where the localized states interact only 
with the electron in the exciton: the full Hamiltonian 
in Eq. is diagonal in the hole spin index /3. This is a 
good approximation for electrons in neutral donors, since 
it is equivalent to neglecting the electron-hole exchange 
interaction which in most semiconductors is much smaller 
than the electron-electron exchange. Concerning the sec- 
ond example we will consider, i.e. the case of the Yb'^+ 
ions in III-V, it is known that these ions act as strong 
isoelectronic traps for electrons and the s-f exchange in 
the conduction band dominates. 

The interaction of the excitons with an external time 
dependent optical field provides the mechanism for the 
control of the two localized spins and is described by the 
Hamiltonian 

HXL - ^ae"'''(bubk=0,a.+f3=a + hc , (3) 

where fl^r is the Rabi energy of the intcrband optical 
transition and tku]^ is the energy of the laser, a is the 
polarization of the light. We have used the rotating wave 
approximation in Eq. (O. 4>is is the envelope function of 
the electron hole pair taken at p = r^—rh = 0. In the case 
of a cw laser field, the time dependence can be eliminated 
by moving to the rotating frame with frequency lol, thus 
replacing ej, by ek — Tiojl in ^o- 

We are deriving an effective Hamiltonian for the two lo- 
calized spins in the presence of the laser field. For a fixed 
value of the polarization a this Hamiltonian is four di- 
mensional, corresponding to the degenerate ground state 
described by |A) = {| TT), I Ti), I it), I ii)} for the two 
localized spins and can be written to the second order 
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perturbation theory in Hxl as 



eff 



1 



XL- 



PHxsP 



Hxl\\') 



(4) 



where P is the projector operator onto the subspace of 
one exciton plus two spins and is the ground state 
energy of two spins with no exciton. Operating Hxl 
onto I A) generates states with one exciton at fc = 
with polarization a, states written as { | XOa) }. The 
ground state energy can be chosen to be =0, and the 
expression above simplifies to 



HxJ = Vj2\^<r\^\ <Pls P Gao.,A'0.(^l) (5) 



where 

GAOo-a'Oo-(t^L) = (AOcr| 



[G"(a>L)]-i-i7i 



lA'Oa) (6) 



is the Green's function operator for the system composed 
by the exciton and two spins, and 



2M . 



We remark that we work in the off-resonance regime for 
which u)L < eoi thus making the real part of always 
negative. Since we are dealing with only single-exciton 
processes, the Lippman-Schwinger equation for G can be 
rewritten in terms of an equation for the T-matrix defined 
by the relation G = Go + GoTGq- We solve the problem 
in two steps: (i) the and operators represent- 
ing the scattering of the exciton with only one impurity 
(identified by the index A or B) is solved, (ii) The T- 
matrix for the exciton interacting with two impurities is 
explicitly rewritten in terms of and usingi& 



T 



1 



1 - T^G^T^G^ 



T^'il + a'T""] + (A ^ B) (8) 



where {A ^ B) stands for repeating the previous term 
with interchange of superscripts A and B. Eq. (jSJ takes 
into account exactly all the multiple scattering processes 
between the exciton and the two localized spins. We 
will focus in the next section on the interaction of the 
exciton with a single localized spin. The role of multiple 
scattering effects in the two spins case are addressed in 
section Hvl 



EXCITON-SINGLE IMPURITY 
SCATTERING 



III. 



This section focuses on the solution of the T-matrix 
equation for one scattering center (named A). Due to the 
short range nature of the exchange interaction, the ex- 
change integral Jk,k' in Eq. (O is often reduced to a con- 
stant, corresponding to a delta-like interaction in space. 



Here we consider a more realistic form of the interaction 
using the separable potential approximationiSiiS where 
Jk.k' = JvkVk', with Vk being a dimensionless form factor 
that depends only on k = \k\. Vk describes the effect of 
the finite size of the non-local exchange interaction. The 
separable form of Jk,k' will allow us to obtain analyti- 
cal expressions for the T-matrix, and provide a flexible 
theoretical framework with parameters that can be taken 
for the experiments. On the other hand, this potential 
can support at most one s-like bound state. The integral 
equation for the T-matrix, T = Hi + HiGqT, can then 
be written explicitly as 

Tkk'aa' yVkVk'S •W + 

+^ E ■ W'GO„r,d,,„„„, . (9) 

kf'a." 

We can write the T-matrix as a sum of a scalar and vector 
part 



kk' aa' 



VkVk' 

V 



-f Tis^ • W] , (10) 



(7) and, using the identity 



3 -s 



(11) 
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we rewrite Eq. (jSJ in the form of two coupled equations 

Ti = J+ JFo(ro-ri/2) (12a) 

To = ^JFoTi (12b) 
16 

where 

^oK)--^E^^"G"-" • (13) 

The reduction of the integral equation to two algebraic 
equations is a consequence of the form of the interaction. 
The two coupled equations in Eq. ((T^ are solved and 
give 



To 
Ti 



3J 







16 1 , 



JFq 

2 
J 



1 



JFq 
2 
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(14a) 
(14b) 



This analytical solution allows us to investigate the 
strong coupling regime in which the quantity JFq is not 
small. The most interesting feature of the strong cou- 
pling regime is the formation of bound states of the ex- 
citon with the impurities, identified by the poles in the 
T-matrix. Varying the frequency of the laser lol, which 
will modify the Fq, we can scan the spectrum to ob- 
tain the energy of those bound states. We remark that 
under the condition lol < co: and assuming reasonably 
that that the potential v{k) is an analytic function of fc, 
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no singularities or branch cuts exist for the function Fq. 
Therefore, the only source of poles is given by the zeros 

of the function 1 + — '^'■'^j^"-' , appearing in Ti and Tq. 
Considering separately the singlet and triplet channels 
we find 



3/4J 



1 + 3/4 JFo 
1 - l/4Ji^o ■ 



(15a) 
(15b) 



In the singlet and triplet channels only one of the two 
poles in Eq. (|14|l is present. Due to the fact that Fo{ujl) 
is negative for all allowed values of wl, we also remark 
that, as expected, the exciton binds in a singlet spin state 
if J > (antiferroniagnetic coupling), while the bound 
state results to be a triplet if the exciton-electron ex- 
change is ferromagnetic (i.e. J < 0). 



IV. EXCITON-TWO IMPURITIES 
SCATTERING 

Starting from the results obtained in the previous sec- 
tions, we construct in this section the solution for the 
exciton-two impurities T-matrix and the corresponding 
Heff for the localized spins. Eq. © can be expanded in 
terms of T operators as 



(16) 



is the corresponding 



The matrix for can be obtained from a simple phase 
shift: if T^f., is the T-matrix for a scattering center with 
potential V{r) then e^*'-'^ ^' 
one for a potential V{r — R) 
matrix elements of Eq. H16|l in the k representation. To 
illustrate how this series can be summed let us consider 
as an example the third term in Eq. H16|l 



■RrpA 

^kk' 

i.e. Tg, We can take the 



(17) 



k",k" 



where we have defined 

Reordering factors and defining the function 
Fn{u^L) = ^J2e^'-^vlGl 



'k J 



The T-matrix is now expressed as an operator in a 8- 
dimensional space generated by three spins ^: one elec- 

(18) tron in the exciton and two localized electron states. By 
direct inversion and products of 8 by 8 matrices, Eq. 
can be rewritten in terms of a combination of spin prod- 
ucts (see Appendix IXjl . and using G — (?o ^" GqTGq we 

(19) obtain the spin dependent effective Hamiltonian 



this term takes the form, 



V 



(20) 



Following the same procedure, the full series can be 
summed to, 



Tk,k' 



VkVk'/V ^'^[^.^^^yS]^^^^-^^) (21) 



1 - FlfT^TB 



(22) 



Bgf f represents an effective magnetic field acting on both 
spins and Jeff is an effective isotropic Heisenberg ex- 
change. The effective magnetic field and exchange con- 
stant can be written as 



and 



eff 



b.sl^vpjl -JF^ ) 



Jeff - 



(1- JF+)[1- JF+(3JF^ -2)]2 



is\^vlJ'Fn/2{l-JFj,) 



(1 - JF+)[1 - JF+{3JFj, - 2)][1 - JFj, {3JF+ ~ 2)] ' 
I 



(23) 



(24) 



where we have defined 



F^{^L)^^Y.^l±e^^-^)vlGl . (25) 



z identifies the direction of propagation of the light, 
(5 = — hiOh is the optical detuning, and f2cr± corre- 
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spond to the contributions to the Rabi energy from the 
two circularly polarized components of the light. From 
Eq. H21|l a spin-independent term is also derived which is 
not shown in Eq. 122|) since it is irrelevant for our pur- 
poses. If we want to include the effect of the degenerate 
light hole band the two expressions in Eqs. (|23|l and H24() 
should be multiplied by 2/3 and 4/3, respectively. By 
keeping the lowest order in J in Eqs. ()23|l and ()24|l we 
obtain 



B, 



eff 



and 



Jeff - 



6^ 



-1-^1 



|2„.2 



vU-+0{J'), (26) 



52 



-\<I>1s\^vP^Fr/2 + 0{J^) (27) 



which recovers the Optical RKKY result of Ref. 0. The 
magnetic field induced by virtual excitons has recently 
been analyzed using a more fundamental approach in the 
case of a single imp urity by Combescot and Betbeder- 
Matibet in Ref. 'M. In this reference the spin indepen- 
dent term that provides a correction to the optical Stark 
shift is also discussed. 



V. SPIN-SPIN COUPLING 

In this section we apply the results obtained above to 
(i) excitons mediating the interaction between two elec- 
tronic spins localized in shallow donors (e.g. GaAs:Si) 
and (ii) excitons mediating the interaction between two 
magnetic ions with spin 1/2 (two rare earth ion Yb'^"'" 
in InP). Yb in InP is one of the most investigated rare 
earth doped III-V material. In principle Yb'^+ in GaAs 
could be used but it is technically more challenging to ob- 
tain samples where only substitutional Yb is present^ 
We will focus on the effect of the binding of excitons 
on the spin-spin coupling. The parameters J and the 
range of the potential Vk can be fixed in such a way that 
the single-spin exciton T-matrix reproduces the binding 
energy and the spin configuration of the bound exciton 
obtained from the experiment. 



A. Shallow donors 

For a scheme of the system we can refer again to Fig.^ 
For excitons interacting with a shallow neutral donor the 
effective mass approximation can be used. The problem 
of excitons bound to neutral donors and acceptors has 
been heavily investigated in the past experimentally and 
theoreticallyi^^ In the case of GaAs it is known that the 
exciton binds to a neutral donor with a binding energy 
of about 1 meV. It is also clear from the magnetic field 
dependence of the bound exciton resonance that the two 
electrons are paired in a singlet around the donor ion and 
the hole is bound by Coulomb interaction. The picture 
is very similar to the one of a positive charge bound to 



an H~ ion. As in the H~ ion case, the dominant term 
responsible for the binding of the two electrons is a ki- 
netic exchange term and we can therefore disregard the 
effect of a spin independent term in the Hamiltonian. 
The range of the kinetic exchange is determined by the 
hybridization between the localized electron state in the 
neutral donor and the electron state in the free exciton. 
We therefore assume that the Vk is of the form 



1 



1 + (AA:)2 



(28) 



where the parameter A determines the range of the po- 
tential. In the following we will use the excitonic atomic 

units where energy is given in excitonic Ry* and lengths 

4 

are in Bohr radius a|j. Ry* — -p^^ where fj, is the re- 
duced mass of the electron hole system (in electronic mass 
units) and Eq is the static dielectric constant in the semi- 



conductor. 



and the relation Ry* — 



2tia*' 



holds. For GaAs these units give IRy* ^ 5 meV and 
a*g ^ 100 A. Using these units we can calculate the func- 
tions Fq and Fr in Eq. and Eq. as 

1 1 



Fo 



Fr 



47rA (AV^ + i^) 



inR 



(29) 
(30) 



where 1/ = fJ-/M is the reduced total mass ratio of the 
excitonic system which is about 1/5 in GaAs taking = 
0.08 and m^, = 0.17. Notice that Fr can be rewritten as 



Fr = Fo Aiy- 



R{Ks/5^ -v) 



(31) 



and has no poles for positive detuning 5\ for i? 3> A, Fr_ 
has a Yukawa form with a detuning-related decay length 
a*Q\fvJb as found in Ref. 6, while at short R the finite 
range of the potential regularizes the 1 jR divergence. In 
Wfc we take A =0.25. Using the fact that the exciton binds 
to the donor only in the singlet channel, we can determine 
the value of J in the in Eq. (|15a|l in such a way to 
have a pole at the experimental binding energy. The J 
is positive, as expected from the fact that the kinetic 
exchange is antiferromagnetic, and we take its value to 
be J = 1 Ry*[a*^Y which gives a binding energy for the 
singlet of 0.23 Ry* ^ in accordance with the experimental 
value of 1 meV. The triplet is unbound. 

We plot in Fig. [21 the coupling constant Jg// obtained 
from Eq. (|24() as a function of the energy of the laser 
measured from the bottom of the free exciton band, 
b — — hujL. A small imaginary contribution to the 
energy, 77 — 0.0001i?y*, has been added in all the plots. 
The Rabi energy is =0.05 meV. The a— component 
of the Rabi energy is zero. The separation between the 
two neutral donors i? is 2 . In the region of large de- 
tuning we have a ferromagnetic coupling in agreement 
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FIG. 2: Coupling constant Je// between the two electronic 
spins localized in a shallow neutral donor embedded in GaAs 
as a function of the laser energy, measured from the bottom of 
the free exciton band. The intensity of the laser corresponds 
to a Rabi energy of Q,a+ = 0.05 meV. The dashed line gives 
the result predicted with the same parameters using second 
order perturbation theory in the coupling constant J. 
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with the results obtained in the ORKKY limit. When 
we approach the energy corresponding to the binding of 
the exciton to the impurity, at 5 = 0.23, we observe that 
the interaction is strongly enhanced and there is a region 
with an antiferromagnetic (AF) coupling. Multiple scat- 
tering between the two impurities results in the forma- 
tion of bonding-antibonding states for the exciton. When 
the light has a frequency in the bonding-antibonding gap 
the effective interaction changes sign. This is analogous 
to the antiferromagnetic coupling generated by superex- 
change in magnetic materials.^'' When the laser is tuned 
above the resonances we recover again the ferromagnetic 
coupling. In the same plot we also show the effective cou- 
pling that would result by keeping the lowest order in J 
(ORKKY). In this case no resonances due to the binding 
of the excitons are present and, in order to obtain a size- 
able coupling, the laser has to be tuned close the bottom 
of the excitonic band. In Fig.|31we show a contour plot 
of the effective spin-spin coupling as a function of the 
detuning and impurity separation. Black corresponds to 
strong negative and white to strong positive coupling, 
the gray tone in the upper-right corner corresponds to 
zero. At large R the coupling is mostly ferromagnetic 
and there is only a small region close to the exciton bind- 
ing energy where the coupling can be AF (white region 
in the plot). When the distance between the two im- 
purities decreases, the bonding-antibonding gap and the 
region corresponding to the antiferromagnetic coupling 
is wider. The dashed line indicates a change of sign of 
Jef f ■ Notice also the different decay of the interaction as 
a function of R for different values of the detuning. At 
S =0A the maximum strength is at R = 0.8 and decays 
quickly within a quarter of a*g to the minimum value in 
the plot. At S =0.1 the same minimum is reached within 
a much larger interval of about 2 a|j. This is consistent 



FIG. 3: Jeff as a function of the donors separation R and 
detuning S . The contour plot identifies the regions where the 
coupling is Ferromagnetic (FM) or antiferromagnetic (AF). 
The (red) dashed lines indicate a change of sign of Jeff- 

with the fact that at a small detuning there is a contribu- 
tion from the free exciton band which can give a longer 
range for the effective interaction. 

B. Rare earth impurities 

The magnetic properties of the Yb'^^ ion in III-\l2^ 
arise from its partially filled 4/ shell, possessing 13 elec- 
trons. In III-V materials, for a substitutional impurity, 
the crystal fields split the ground manifold of the ion 
into two doublets (spin=i), Fg and F7, and a quadru- 
plet (spin=|), Fg. The lowest lying state is the Kramers 
doublet Fg, which behaves like a spin ^ with an effective 
isotropic g = 24/7i^ Yb in InP replaces Indium and acts 
as an isoelectronic trap. From electrical^ and optical2& 
measurements it is known that the exciton binds to this 
isoelectronic impurity with a binding energy of 30 meV. 
The binding is due to a short range potential that arises 
from the difference in the core pseudopotential between 
the impurity and the host ion it replacesj22i It is reason- 
able to assume that this short range potential is spin- 
independent and we take it into account by adding to 
the exciton-impurity Hamiltonian of Eq. the term 

i/2 = -^E^M-'(l + e-^('='-'^-)^)6L„^5fe„^ . (32) 

k,k' 

We use the separable potential approximation also 
for the spin-independent short-range potential and we 
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parametrize it in the form A^ fc' = AiifeV^/ , i.e. it has the 
same k dependence of Jk,k'- A more general analytical 
result can be obtained using a separable form for Ak^k' 
with different coefficients, but we expect the range of the 
s-f exchange and that of the impurity potential to be 
very similar. The value of A is determined by imposing 
that the exciton-single impurity T-matrix has a pole for 
both singlet and triplet channels at 30 meV. Following 
the same procedure used in Sec. Illll we obtain for the T 
operators in the singlet and triplet channels 



rpS 



-3/4J + A 



1 + 3/4 JFo - AFo 
^ J/4 + A 

1 - f/4JFo - AFo 



(33a) 
(33b) 



The expressions for the Je// and Be// modified by the 
presence of A can be obtained by plugging the Eqs. 
in the general expressions of Eqs. (|A3ll in Appendix ^ 
The quantity J is the s-f exchange interaction between 
the impurity and the electron in the exciton. In typi- 
cal rare earth ferromagnetic semiconductors the s-f ex- 
change is ferromagnetic and is of the order of few eV 
A'^ comparable to the s-d exchange in Mn based 
diluted magnetic semiconductors 3^ We are using J = 
— 10"'' in our units which corresponds to a conservative 
estimate of 0.7 eV in InP. In InP the value of the Ry* 
is about the same as that of GaAs (5 meV), while the 
Bohr radius is about 120 A. For A in we take A = 0.01 
which is of the order of the ionic radius of Yb'^^ . 

We show in Fig.^the contour plot of Je/ / as a function 
of the laser detuning S = eo — ^ojl and of the separation 
between the impurities R. At large distances we observe 
two resonances related to the binding of the exciton in the 
singlet and triplet channels. Fig.|Sl(a) shows in detail the 
Jeff for a distance R= 1 a^. The two peaks in Fig.[Sl(a) 
correspond to the exciton bound to the impurity in the 
triplet and singlet channel. The peak at larger detun- 
ing corresponds to the triplet since the s-f exchange is 
ferromagnetic. For shorter distances we see from Fig. 0| 
that each of the two peaks starts to split. The singlet 
(at smaller detuning) follows a behavior similar to the 
one of the shallow donors described above: the bonding 
and antibondig states identify a region where the cou- 
pling become antiferromagnetic. The triplet state splits 
in many different peaks as can be seen from Fig. [S] (b). 
The sign of the interaction can change many times as 
a function of the detuning in this short distance region. 
This is indicated by the sign of Je/ / plotted in the lower 
part of Fig. |5l (b). Overall the antiferromagnetic cou- 
pling dominates at short distances while the interaction 
is ferromagnetic at large R. 



VI. DISCUSSION 

The spin-spin interaction control discussed in this pa- 
per has potential applications in quantum computing im- 
plementations. In fact, an optical control of the spin of 
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FIG. 4: Coupling constant Je// between the two magnetic 
Yb"^"'" localized in InP as a function of the laser detuning 
5 = e — TujJl and separation between the ions. The dashed 
(red) lines indicate J^ff—O and a change form ferromagnetic 
to antiferromagnetic coupling. The intensity of the laser cor- 
responds to a Rabi energy of SI = 0.1 meV. 



electrons localized in quantum dots or impurities has sev- 
eral advantages with respect to approaches where elec- 
trodes are needed. Ultrafast lasers are available, promis- 
ing the realization of quantum gates in time scales that 
are hard to achieve with an electrical control. Laser are 
also very flexible from a control perspective since pulse 
shaping can be used to advantage in the accuracy and 
speed of quantum operations »^ Finally, metallic elec- 
trodes necessarily add a source of noise for the quantum 
system, and they are not needed in an optical scheme. 
The possibility of changing the sign of the spin-spin inter- 
action exploiting the exciton bound states can add flex- 
ibility to many control schemes for the qubits, like e.g. 
in the only exchange scheme."^' We have seen that reso- 
nances in the spin-spin coupling induced by the binding 
of the excitons can increase the magnitude of the inter- 
action for distances that are reasonable from a nanofab- 
rication point of view. This will imply that lasers with 
lower intensities can be employed in the control. The 
polarization of the light represents an additional control 
parameter that can be used to selectively address qubits 
with an optically induced magnetic field. This is also an 
advantage from a practical point of view since it could 
simplify the experimental setup by eliminating the need 
of an external magnetic field. 

Although the feasibility of single impurity spectroscopy 
in semiconductors has been proven'^^''^^, little attention 
has been paid to optical properties of impurity-bound 
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FIG. 5: Coupling constant Je// between two Yb ions in InP 
as a function of the detuning (a) Large distance. The cou- 
pling is ferromagnetic and the resonances in the interaction 
are close to the energy of the exciton bound to the Yb. (b) 
Short distance. The triplet channel splits in many different 
peaks producing many changes of the spin-spin coupling sign. 
The lower curve shows the sign of the coupling constant. 



excitons for information storage and processing. Impuri- 
ties deserve at least the same attention as quantum dots 
for such applications. Their homogeneous character and 
the variety of properties that one can obtain combining 
different hosts and ions are indeed pecuHar advantages. 
An exciton bound to an impurity has optical properties 
very similar to an exciton trapped in a shallow quantum 
dot. Most of the ideas involving excitons in quantum 
dots as a main ingredient for quantum information and 
communication can be reformulated for excitons bound 
to impurities. We have provided only two examples here, 
but our phenomenological theory, being based on inputs 
from the experiments, is very flexible and many other 
combinations of host and ions can be used to explore a 
large range of confinement energy and different optical 
properties. We also have seen that the spin-spin cou- 
pling has a resonant behavior at frequencies depending 



on the separation between the impurities. By organizing 
the impurities in chains with different separation this can 
be used to selectively address a single pair of impurities 
and it allows for scalability. 

A very special case is represented by impurities in sili- 
con. This material has obvious technological advantages 
and many proposals for using impurities in Si for quan- 
tum computing have been suggested i^^iSiS^ In particular, 
the optical control of electronic spins localized by deep 
donors in Si using a control impurity has been proposedi^ 
In the scheme we are suggesting here, the exciton bound 
to the impurity plays the role of the control impurity 
and it takes advantage of the host material for mediating 
the interaction. Even if Si is an indirect gap material, 
there is a finite optical coupling to the exciton bound 
to the impurity due to symmetry breaking. Additional 
complications in the use of excitons bound to donors for 
mediating spin-spin coupling arise form the valley degen- 
eracy in Si."^^ We will address the optical spin control of 
impurities in Si and the role of valley degeneracy in a 
future publication. 

Excitons bound to rare earth magnetic ions can be 
controlled very rapidly and efficiently due to their strong 
dipole moment. Their dipole moment is mainly deter- 
mined by the optical properties of the host material, 
since it involves the creation of electron hole pairs across 
the semiconductor gap. At the same time, they inter- 
act with the internal degrees of freedom in the core / 
states. Schemes involving excitons bound to rare earth 
impurities in III-V materials bring in the advantages of 
the optical properties of the host and the stability of the 
internal degrees of freedom of the / orbitals in the rare 
earth ion where the qubit is stored. This hybrid system is 
thus extremely powerful, providing both reliable storage 
and fast processing of information. 

Finally, the light controlled spin-spin coupling in a 
semiconductor matrix is also appealing for the coherent 
control of macroscopic properties of materials. This was 
the idea behind the coherently induced ferromagnetism 
m There, a finite critical temperature for a para- 

magnetic to ferromagnetic transition in diluted magnetic 
semiconductors was found when the material is coupled 
to a strong laser field. The results presented in this paper 
suggest that the presence of bound states could enhance 
the effect. Also, the same idea could be used in other 
systems where the light can induce antiferromagnetic or 
glassy phases starting form a paramagnetic system. This 
represents a unique opportunity to study phase transi- 
tions in a solid where the coupling is controlled by an 
external field and may open to a new class of controlled 
materials to be investigated. 

In conclusion, we have studied the problem of two 
spins ^ localized by impurities in semiconductor in the 
presence of an intense light field. The light induces a 
frequency dependent spin-spin coupling and a magnetic 
field that can be controlled by the polarization of the 
light. The effects are enhanced by the presence of im- 
purity bound excitons which may split in bonding anti- 
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bonding states in the case of two impurities. The sign of 
the spin-spin couphng is generaUy ferromagnetic, but it 
can switch to antiferromagnetic when the laser is tuned 
in the bonding-antibondig gap. We have developed a 
flexible theoretical approach based on scattering theory 
where the parameters from the experiment can be used 
to estimate the size of the effect. We have discussed ex- 
plicitly the case of two neutral donors in GaAs and two 
rare earth magnetic ions (Yb'^+) in InP. 



-PrT" ^] + {A ^ B) . The traceless part of this matrix is 
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APPENDIX A: MATRIX REPRESENTATION OF 
THE r OPERATOR 



where a and /3 can be conveniently expressed as a func- 
tion of the single impurity and operators in 
Eqs. or Eqs. ((23) as 



/3 = 



{T^ ~ T^)Fr 



2FbT'^TS) - 2 



;A3a) 



^ {Fr{TT -TS + 2FrTTTS) - 2) 



(A3b) 



where we have dropped the {A,B} index since we are 
considering two identical centers. Notice that the matrix 
in Eq. (|A2ll can be rewritten as 



Using the basis set |s^, sf , s^) we obtain for the spin 
products ■ s, ■ s, ■ the matrices 
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By substituting these expressions in T"^ and and then 
in Eq. H21(l , we obtain after matrix inversions and multi- 
plications an expression for T = (1 — F^T'^T^)^^T"^[1-|- 
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